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We present a time-dependent perturbative approach adapted to the treatment of intense pulsed 
interactions. We show there is a freedom in choosing secular terms and use it to optimize the 
accuracy of the approximation. We apply this formulation to a unitary superconvergent technique 
and improve the accuracy by several orders of magnitude with respect to the Magnus expansion. 
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Perturbation theory when combined with a specific 
treatment for resonances is quite well understood in clas- 
sical and quantum mechanics for time-independent sys- 
tems. This includes also time-periodic driven systems for 
which the periodicity can be treated by Floquet theory in 
a way that yields a time-independent formulation [l|, |3 ■ 
One knows that resonances yield divergent terms, that 
appear as small denominators, which have to be specif- 
ically removed. The counterpart of the concept of reso- 
nance for time-dependent systems is generally associated 
to secular terms whose size grows with time (see Ref. Q 
and references therein). 

With the advent of short (~ 10 fs) and intense (10^'^~ 
10^^ W/cm^) laser pulses, atomic or molecular systems 
can be strongly perturbed in a timescale shorter than 
characteristic times corresponding to the free evolution 
of the system and adiabatic theories are not applicable 
(see, e.g., Ref. |^). The goal of this paper is to formulate 
a time-dependent perturbation theory well adapted for 
perturbations localized in time. 

The conceptual framework of perturbation theory can 
be described as follows: The Hamiltonian of the consid- 
ered system can be decomposed as the sum of two terms 
Hi = Hq + eVi . The first term iJo is assumed to have a 
structure simple enough to lead to explicitly known solu- 
tions for its associated propagator UHo{t,to). The term 
eVi is supposed to be small with respect to Hq, in a sense 
specified below. Time-independent perturbation theories 
can be equivalently formulated at the level of eigenvec- 
tors or operators (5|. A large class of these approaches 
amounts to construct a unitary transformation T such 



that 



T^HiT = H'' + e'V', 



(1) 



where H° is still of simple structure [i.e., its propagator 
UH^it^to) can be explicitly computed] and e'V is a per- 
turbation whose size is smaller than the original one. To 
compute the transformation T explicitly, one represents 
it in general either (i) in terms of some power series 



(2) 



or (ii) by an iterative construction as a composition of 
transformations 



T 



(3) 
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These procedures generally differ. The former one is re- 
ferred to as the time-independent Poincare-Von Zeipel 
technique, which has been shown to be equivalent to the 
usual Rayleigh-Schrodinger perturbation theory [^. The 
latter procedure includes the Van Vleck technique (for 
which Efc = e'^) and the superconvergent Kolmogorov- 
Arnold-Moser (KAM) expansion (where Cfc = and 
W^fe is Efe-dependent) Q. The perturbative procedure con- 
verges if the remaining perturbation e'V can be made to 
go to zero, as the number of terms in the power series 
^ or as the number of compositions in Eq. Q goes to 
infinity. 

In this description one has to state precisely what class 
of Hamiltonians H'^ can be considered simple. For the 
first order or the first iteration, one considers H'^ ~ Hq + 
eDi with the condition that Di should be compatible with 
Hq in the sense that if the propagator of Hq is known, 
that of Hq + eDi can also be obtained explicitly. In the 
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case of time-independent Hamiltonians the condition of 
compatibility is 

[Ho,Di] = 0. (4) 

For the case of time-dependent Hamiltonians, we show 
that the condition of compatibility can be generalized to 

[H,{t),D,{t)]=^^- (5) 

The construction of transformations of the type of Eq. 

or involves finding the generator —i€Wi{t) of the 
transformation Ti {t) , that is the solution to 

I [W^{t),H^{t)] + V^{t) -D,{t)^^. (6) 

This equation, together with the constraint Q, are usu- 
ally called cohomology equations in the time-independent 
case and are here generalized to the time-dependent 
case. These cohomology equations have the same form 
for higher orders or successive iterations. 

Here we formulate the time-dependent perturbation 
theory by transforming directly the evolution operator 
instead of considering the perturbed Hamiltonian as is 
usually done in time-independent theory. We obtain per- 
turbative corrections to the full propagator in the form 
of a product of propagators which exhibit free parame- 
ters appearing through the general solutions of related 
differential equations. We recover in particular the Mag- 
nus expansion Q as a special case of the time-dependent 
Poincare-Von Zeipel theory. This extension also gives the 
precise correspondence between time-independent reso- 
nances and time-dependent secular terms. In the context 
of pulsed perturbations with a finite duration, the secular 
terms need not be eliminated. We show the remarkable 
result that they can be used to improve the convergence 
of the method at a given order. This optimization is 
achieved without any a priori knowledge of the solution 
by locating the minimum of a given eigenvalue as a func- 
tion of the relevant free parameters that are identified. 
The efficiency of the method is illustrated on a two-level 
system driven by a short intense pulse. 

Perturbation theory, resonances and secular terms. We 
consider the Hamiltonian Hi{t) — Ho{t) + eVi{t), where 
Ho{t) is associated to a known propagator UHo{t,to). 
The formulation is presented here for the KAM method, 
consisting in iterations of transformations which are ex- 
actly of the same form at each step. The first iteration 
involves a unitary operator Ti(t) which transforms the 
propagator UHi{t,to) according to 

Tl{t)UHAtM)Ti{to) = UnAiM, (7) 

into a propagator [///^(^i^o) associated with the sum 
H2{t) of an effective Hamiltonian Hf{t) = HQ{t) + eDi{t) 
which contains contributions up to order e and a remain- 
der e^V2(i). This new propagator, generated by a sum of 
two Hamiltonians, can be written as the product 

UH,{t,to)^UHl{t,tQ)R2{t,to), (8) 



where R2{t,to) satisfies the Schrodinger equation whose 
Hamiltonian is e^UH°{to,t)V2{t)UH°{t,to)- Similarly 
we can factorize UHi{t,to) — UHo{t,to)Si{t,to), where 
Si{t, to) is a unitary operator related to eDi{t): 

i^^Siit,to) ^ eUH„ito,t)Di{t)UHoit,to) Si{t,to). (9) 
The full propagator reads 

UHAt,to)^Ti{t)UHo{t,to)Si{t,to)R2{t,to)Tl{to), (10) 

which yields the first order KAM approximation for 
R2{t,to) replaced by the identity. In this construction 
the only restriction on the self-adjoint operator eDi{t) is 
that it be of order e. Hence we have the freedom to choose 
the Hamiltonian in Eq. as ^-independent, giving 

Di{t) - UH,{t.to)Di{ta)UHo{to.t), (11) 

with Di{to) arbitrary, which is the general solution of 
Eq. 0. This allows one to obtain the solution of Eq. 
® as 

5i(t,to) -exp[-i(t-io)ei^i(to)]. (12) 

Differentiating Eq. H1U|) and substituting Ti{t) = 
exp (— zeWi(t)) leads to Eq. ©. The rest involves a 
series of k nested commutators that reads e^V2{t) = 
EZi JkTW''''^' ' • • • it),kV^it)+D,it)]■■■]. 

It has exactly the same structure at each iteration which 
is useful for applications, particularly when high-order 
computations are needed. 

Iterating the time-dependent KAM algorithm reduces 
the size of the remaining perturbation in a superconver- 
gent way from order to at step n. The time- 

dependent Van Vleck technique would allow one to re- 
duce the size of the remainder from order to e'^'^^. 
These methods, in the formulation presented here, are 
unitary upon truncation. The superconvergent character 
of the KAM algorithm has been shown numerically by 
applying the method to a two-level system perturbed by 
a short time-dependent interaction |lO| . 

The time-independent problem., i.e., the problem of 
finding a transformation Ti that enables one to sim- 
plify the time-independent Hamiltonian Hi according to 
T^HTi = Hq + eDi + e^V2, is recovered when one conve- 
niently chooses Ti as time- independent. In this case all 
the operators, and in particular Di and Wi, are time- 
independent and the standard cohomology equations are 
recovered: [Ho,Di] = Q and Vi - Di + i [Wi,Ho\ = 0. 
Their solutions can be determined using the following key 
property 0]: Wi exists if and only if IIho{Di — Vi) = 0, 
where T^Hq is the projector in the kernel of the applica- 
tion A i-^ [A, Hq] (for an operator A acting on the same 
Hilbert space as H^). The projector Hhq applied on an 
operator A captures thus all the part B of A which com- 
mutes with Hq: [B^Hq] = 0. The unique solution Di 
allowing Wi to exist and satisfying Eq. Q is thus 

Di^Uh,Vi= lim i r e-''"-Vie''"\ (13) 
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The resonances are associated with terms of Vi which 
commute with Hq. Apphcation of Eq. can be inter- 
preted as an averaging of Vi with respect to Hq which 
ahows one to extract resonances. 

For the time- dependent problem, the general solution of 
Eq. © reads (up to a term UHo{t7to)BiUHo{tQ,t) with 
Bi any self-adjoint operator that we set here to 0) 

Wiit)^ [ dsUH„it,s){Vi{s)-Di{s))UHois,t), (14) 

with t[ any real number. Defining the average 
1 /■* 

n_Fi= lim - / dsUHo{t,s)Vi{s)UH„{s,t), (15) 

one can show the following property: if Wi{t) is bounded 
for negative infinite times, then n_(Vi — Di) = 0. This 
is satisfied by Di = II^Vi, the only solution compatible 
with Eqs. (0) and 1)15(1 . Hence, the averaging _Di ~ H-Vi 
allows one to remove secular terms at negative infinite 
times. This gives the precise correspondence between 
the resonances of stationary problems and the secular 
terms of time-independent problems. We remark that 
the definition H15|) of the average can be in fact recovered 
from the formal calculation of the average n/f^T^ [cf. 
Eq. H13|l ] with respect to Kq — —i-§f + Hq in an extended 
space, which includes time as a coordinate [lol[llj|. 

For a pulsed perturbation that is switched on at 
ti and off at tf, Eq. becomes n„Vi = 

UHo{'t,ti)y{ti)UHo{ti,t)- This is a particular so- 
lution of Eq. © corresponding to the choice 
Di{to) = UHoito,ti)Vit^)UHo{U,to) in Eq. (EJ. 
An alternate definition of the average: n^.^! = 
lim^^oo i /t*"^^ dsUffg {t, s)Vi {s)Uh„ (s, t) gives a different 
averaging II+Vi — UHo{t,tf)V{tf)UHa{^f,t) and allows 
one to remove secular terms at positive infinite times. 

Generally one cannot remove simultaneously the sec- 
ular terms at negative and positive large times. This 
shows a conceptual difference between stationary reso- 
nances and secular terms associated with perturbations 
localized in time. Furthermore, it appears that the av- 
eraging such as Eq. I|15|) is not appropriate, but that a 
definition which combines the two definitions gives a new 
secular term that could improve the convergence of the 
algorithm. This suggests to work with the general so- 
lution Hll|) of Eq. (O, written with the perturbation 
evaluated at a free time ti as the arbitrary operator: 

Di{t) = UHo{tM)Vl{tx)UHo{tl.t). (16) 

The free ti can then be chosen to minimize the rest af- 
ter the first iteration, as we describe below. One has n 
such free parameters i^, k = l,rt for n iterations of the 
KAM algorithm. There is only one such free parame- 
ter for the time-dependent Poincare-Von Zeipel and Van 
Vleck methods that are order by order techniques. An 
interesting result is that we recover the Magnus expan- 
sion from the time-dependent Poincare-Von Zeipel in the 
particular case of Dk = and t'f. = to for k = l,n. 



Optimization of the perturbation theory. After one it- 
eration, the rest i?2(i,io) defined in Eq. |(SJ) is associ- 
ated with a second order operator through i?2(i,^o) = 
g-ie^G2(t) ^[^i^ G2{to) = 0. The closer i?2(<,io) is to the 
identity, the smaller the correction terms are, i.e., the 
more accurate the approximation is. We evaluate the 
lowest order contribution to e^G2{t) as 

e^G^^\t) = £2 /■ duUHiito,u)V2{u)UHfiu,to). (17) 

Jto 

It is this operator that has to remain small for the al- 
gorithm to converge. The size of an operator A can be 
characterized by the norm ||A|| = sup||^||^]^ ll^'0ll with 
■0 in the appropriate Hilbert space. For an Hermitian 
matrix this norm reduces to the largest of the absolute 
values of its eigenvalues. 

In order to improve the accuracy we thus seek to mini- 
mize A2 (t) , the largest of the absolute values of the eigen- 
values of e^G'2^''(i), with respect to the free parameters. 
To optimize the KAM algorithm, we have at our dis- 
posal two free parameters tk and tj, at each iteration. 
We expect the parameters tk to significantly affect the 
convergence, as they are related to secular terms. 

Perturbation theory for short intense pulses. We con- 
sider a system described by the Hamiltonian H (au- 
tonomous or not) and perturbed by a time-dependent 
Hamiltonian V{s) whose characteristic duration is r. The 
perturbation is assumed to satisfy [l^(s), t^(so)] — 0, 
Vs, So which is realized in many situations of interest. 
We define a sudden parameter e as follows. A dimen- 
sionless time t and dimensionless operators H and V{t) 
are defined through s = rt, H = huiH, and V{s) = 
^V{t), leading to the dimensionless Schrodinger equation 
i-§-^U{t, to) = {V{t) + eH} U{t, to), where the sudden pa- 
rameter is defined as e = ut. We then apply the per- 
turbation theory described above with the identification 
Ho{t) = V{t) and Vi = H. This formulation is suited to 
treat intense short pulses. 

Illustration on a pulsed-driven two-level system. 

We consider the case where Ho{t) = n{t)ai and Vi = 
(73 with Q,(t) a pulse that is switched on at ti and off at 
tf, and (Tfc the Pauli matrices. Notice that, as discussed 
above, the role of the perturbation and reference Hamil- 
tonian is interchanged. The pulse area A = /;/ 0(u) du 
is a dimensionless parameter that can be fixed indepen- 
dently of the sudden parameter e. The error between the 
numerical solution of the Schrodinger equation at the end 
of the pulse and the result of n iterations is defined as 
A„ = II Uhi {tf,ti) — u'^^{tf,ti)\\. We use the pulse shape 
n{t) = 2Asin^ {-Kt) for < t < 1, and elsewhere. 

The upper panel of Figure 1 shows a comparison of the 
error Ai for the first order Dyson, first order Magnus and 
one-iteration KAM methods as a function of ti , for a non- 
perturbative area chosen to produce comparable errors 
Ai for the Magnus and non-optimized (ti = 0) KAM 
techniques. The lower panel displays A2 the largest of the 
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FIG. 1: (a) Common logarithm of the error Ai for the first 
order Dyson expansion (dashed line), the first order Magnus 
expansion (dotted line) and the first KAM iteration (solid 
line); and (b) largest eigenvalue of e^G2^'(tf) as a function of 
ti, for A = 1, e = 0.5 and t[ = 0. 

absolute values of the eigenvalues of e'^G^\t[) defined in 
Eq. (fTTji . We clearly see that the error of the first KAM 
iteration is correctly estimated by this eigenvalue A2 and, 
in particular, minimized when A2 is minimized, i.e., for 
the value t\. It is worth noting that modifying ti covers 
more than one order of magnitude in the error, a situation 
that is not restricted to these values of the parameters. 
The optimized solution provides an improvement of the 
accuracy by almost one order of magnitude with respect 
to the Magnus calculation. 

Figure 2 displays a comparison of the error A2 for 
the second order Dyson, second order Magnus and two- 
iteration KAM methods as a function of t2, for ti — t*. 
It is seen that the Dyson approach is not applicable in 
this context of strong field as the second order performs 



worse than the first one. Figure 2 also shows that the 
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FIG. 2: Common logarithm of the error A2 for the second or- 
der Dyson expansion (dashed line), the second order Magnus 
expansion (dotted line) and the second KAM iteration (solid 
line) with the same parameters as Fig. 1 and ti — t* ~ 0.39. 

second KAM iteration can be enhanced by about two or- 
ders of magnitude with an appropriate choice of ti and 
t2- This optimized second KAM iteration provides an 
improvement by two and a half orders of magnitude with 
respect to the second order Magnus technique. Higher 
iterations of the KAM technique can also be optimized 
and produce still better improvement owing to its super- 
convergent character. 

In conclusion, we have presented an optimized pertur- 
bation theory for pulse-driven systems, which applies to a 
wide class of processes controlled by intense femtosecond 
laser pulses. The optimization reduces to the evaluation 
of eigenvalues and is therefore easy to implement. We an- 
ticipate that this approach will be useful! in the context 
of the laser control of atomic and molecular processes, 
such as the phase space localisation of Rydberg electron 
[T^. or the alignment and orientation of molecules 0|. 
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Photonique and Conseil Regional de Bourgogne. 
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